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Theory of the Effect of Drag on the Orbital 
Inclination of an Earth Satellite 

John P. Vinti 

The rotation of the earth's atmosphere with the earth gives rise to a change in the 
inclination of the orbit of an earth satellite. In this paper it is assumed that the drag is in 
the direction of the air velocity relative to the satellite and that its magnitude diminishes 
so rapidly with altitude that it is appreciable only near perigee in an elliptic orbit. With- 
out further assumptions about the properties of the drag function, results are then deduced 
for the secular changes in inclination up to, but not including, the final ballistic stage of 
motion. 

1. Introduction 

Although the drag on an earth satellite is not well known, it is possible to express the 
secular changes that it produces in the inclination of the orbit, in terms of easily measurable 
quantities, provided only that the drag is given by a certain very general law. The present 
paper, which was suggested by a note by Wild hack, 1 gives a more complete treatment of the 
problem, with more liberal use of mathematical methods. 

2. Forces and Frame of Reference 

If we adopt a Cartesian reference system with origin at the center of the earth and axis 
directions fixed relative to the fixed stars, the following forces act on a satellite: 

F e =inverse square gravitational force of the earth, 
F = correction to F arising from oblateness of the earth, 
F s = gravitational force of the sun, 
F m = gravitational force of the moon ; 

F f5 =inertial force arising from acceleration of the earth, produced by the sun, 
F im =inertial force arising from acceleration of the earth, produced by the moon, 
F rf =drag force. 
There are other forces, all of which may be disregarded as extremely small. Moreover, 
F, v +F t - S and F TO +Ff TO are small compared to F . We are thus left with the forces F,., F , and 
F,. 

The oblateness force F has a secular effect on the orbit, inasmuch as it produces a secular 
precession of the line of nodes and of the line of apsides. In considering changes in the in- 
clination of the orbit, however, we are concerned with nutation and the effect of oblateness is 
only to produce a small oscillatory nutation. We shall therefore be concerned only with the 
forces F e and F d , up to that point where we find results dependent on the (oblateness-induced) 
motion of the perigee. At that point it will then be sufficient to quote results by Garfinkel 2 
to find how oblateness influences the secular effect of drag on the inclination. 

3. Fundamental Theory 

With the above choice of reference system the center of the earth plays the role of a fixed 
point and the angular momentum p of the satellite is given b}^ 

p=m rXr. (1) 

Here m is the mass of the satellite, r is the radius vector to it from the center of the earth, 
and r is the velocity of the satellite. 

i W. A. Wildhack, Science, 128, 309 (Aug. 8, 1958). 
2 B. Garfinkel, Astron. J. 63, 88 (March 1958). 
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As an orthonormal vector basis we choose k pointing northward along the axis of rotation 
of the earth, i pointing in such a fixed direction that at time t=0 its extension pierces the 
Greenwich meridian, and j so that j = kXi. 

The orbit is a space curve, but it is natural to define the instantaneous plane of the orbit 
as the plane of two successive radius vectors, i. e., as the plane of r and r or the plane perpen- 
dicular to p. It is then convenient to define the directed unit normal to the orbit as a unit 
vector l p along p. The sense of the motion of the satellite relative to l p is then given by the 
right-hand screw rule. We then define the inclination of the orbit as the angle between k 
and l p . (In terms of observable quantities is thus the greatest latitude reached by the sat- 
ellite in a revolution.) Then 

V=plv, (2) 

where p= jp| and 

l p =i sin cos </>+j sin sin <£+k cos (3) 

the angle </> being the angle between i and the projection of p on the equatorial plane {i, j}. 
Then p, 0, and <£ are the polar coordinates of p. 

The vector 

l^kXl? csc (4) 

is a unit vector lying both in the equatorial plane and in the plane of the orbit, so that it points 
along their intersection, the line of nodes. If l r is a unit vector along r, it follows from (1) 
that l r is perpendicular to \ v . Moreover 1^ is perpendicular to 1^, by (4), and so is \ v YAn> 
Thus the unit vectors l r , 1^, and l p Xljv are coplanar, since they are all perpendicular to l p . 
If now yp is the angle between the radius vector and the line of nodes, so that 

*=(1*, lr), W 

we have 

l r =liv cos ^+lpXljv sin ^, (6) 

where the positive sense for increase of \[/ is from 1^ towards l^Xljv- 
From (3) and (4) 

1^=— i sin </>+-j cos <j> (7) . 

kXl;v=— i cos </>— j sin <£. (8) 

Comparison of (8) with (3) shows that 

l p =^k cos 0— kXl.v sin (9) 

By (4) it follows that 1^ is perpendicular to k, so that k, 1 N , and kXliv form another right- 
handed orthonormal basis, which turns out to be useful. 

4. Changes of Angular Momentum 

If F is the total force acting on the satellite, then 

p=rXF. (10) 

The secular changes with which we are here concerned are produced only by F d , so that for 
our present purposes 

P=rXF,. (11) 

To find the rates of change of p, 0, and </>, it is now convenient to express each side of (1.1) in 
terms of the basis k, l N , and kXW- 
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The left side may be so expressed by menus of (2), (3), (7), (8), and (9), with the result 

p=k(p eos d—p sin 6)6+l N p sin 0<£+kXlv(— $ sin 6 — p cos 00). (12) 

5. Drag Force and Drag Torque 

We assume the atmosphere of the earth to !>e rotating with il like a rigid body. If to, is the 
angular speed of rotation of the earth, the velocity w of the atmosphere (w for "wind")* relative 
to the inertial system i, j, k, is 

w=o> e kXr. (13) 

The velocity of the satellite relative to the inertia! system is r. Its velocity v a relative to 
the atmosphere is then 

v a =r— w. (14) 

We assume the drag force to be in the direction of — v„, with magnitude dependent on r and 
the relative speed jr— w|. Thus 

F d =/(r,|f-w|)(w-r). (15) 

Equation (15) is the "very general law" mentioned in the Introduction. If the drag* is quad- 
ratic in the air speed |f — w|, then the function/ is proportional to the air speed. No such 
special assumption is needed, however, for the purposes of this paper. 
The drag torque is now given by 

rXF d =/rXw-/rXr. (16) 

Now by (1) 

rXr=p/m, (17) 

by (13) 

rXw=co e r 2 [k-l r (l r .k)], (18) 

and by (6) and (9) 

l r = k sin sin ip+l N cos t/' + kXlv cos sin \p. (19) 

From (16) through (19) it then follows that 

rXF d = - /p/m+/av< 2 lk(l-sin 2 sin 2 yf/) — l N sin sin i£ cos ^— kXl* sin cos sin 2 fl. (20) 

'Phis is the expression for the drag; torque. 

6. Basic Equations 

On inserting (12) and (20) into (10) and equating corresponding coefficients of the basis 
vectors k, l v , and kXlv, we find 

p cos 0— p sin dd=—fm~ l p cos 0+/co e r 2 (l — sin 2 6 sin 2 \p) (21.1) 

p sin 6+p cos 66=—fm~ 1 p sin 0+/o?,r 2 sin cos 6 sinV (21.2) 

p sin 0<£=— ju e r 2 sin sin \j/ cos \p. (21.3) 
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Multiplication of (21.2) by i=^—l and addition of the resulting equation to (21.1) gives, after 
division by expid: 

p J rip6=fm~ 1 p-{-JQi e r 2 [cos 2 \f/ exp(— id)-}- cos 9 sin 2 \//], (22) 

from which there follow 

p=frn- 1 p+fu e r 2 cos (23.1) 

p6= — fu e r 2 cos 2 ^ sin 6. (23.2) 

Equations (23), together with (21.3), are the basic equations for the changes in p, 6, and <t> 
that arise from the drag. 

7. Stages of the Motion 

We assume that the motion of a satellite takes place in three stages: 

Stage (1) : The motion is that of a precessing and shrinking ellipse, with diminishing 
eccentricity, in which the elliptic elements are slowly varying functions of time. For our 
present purposes the relevant elliptic elements I, e, and \[/ p , along with the related quantities 
p, r p , and r a , may be chosen as follows: 

p= magnitude of the orbital angular momentum 
Z=semi-latus rectum 
e= eccentricity 
r p =radius at perigee 
r a = radius at apogee 
\/y p = value of ^ at perigee. 
The elliptic orbit is then characterized by the equation 

r=l[l + e cos^-^)]- 1 , (24) 

where 

l=(GMm 2 )-y (25) 

and 

p = mr 2 \j/, (26) 

M being the mass of the earth. From (24) it follows that 

r v =l{\ + e)~' (27) 

r a = l(l-e)- 1 (28) 

from which 

l=2r v r a (r a + r p )- 1 (29.1) 

e=(r a -r p ){r a +r p )-K (29.2) 

Stage (2) : A quasi-steady spiral stage, in which there are no apsides, so that r is small 
and negative throughout the stage, and in which the speed v is given approximately by the 
formula for a circular orbit without drag. With neglect of oblateness this relation is 

v 2 =GM/r. (30) 
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Stage (3): A final "ballistic" stage, which consists of a rapidly diminishing spiral, and 
in which v is not given approximately by the formula (30) for a circular orbit. Quantitative 
treatment of such a final stage would require an accurate knowledge of the drag function, so 
that we omit consideration of it in this paper, restricting our attention to the changes in 

inclination that may occur in the first two stages. 

We let r P oandr a o be the initial perigee and apogee radii in stage (1) and/', their final common 
value at the end of stage (1). The quantity r x is then the initial radius in stage (2), which we 
take to have a final radius r 2 . In any actual case the quantities r p0 and r a0 would be given as 
initial conditions, while the radii f\ and r 2 would have to be determined by observation — or 
calculated with use of an accurate expression for the drag. For rough estimates of the order 
of magnitude of the expected changes in inclination we shall here simply postulate arbitrary 
values of r p0 and r aQ and try to guess plausible values of r x and r 2 . 

8. A Simplification 

Oblateness is known to produce no secular change in the inclination 0, but only an oscilla- 
tion with amplitude less than 0.05°. Suppose we also anticipate a result of the present paper, 
viz, that drag produces only a small secular change in 6, of the order 0.3°. For a polar orbit, 
with O =9O°, we can then neglect the term in (23.1) that involves cos 6. 

For any other orbit this term, although nonvanishing, is still small in comparison with 
the term/m -1 y;. To show this, note that the ratio R of the two terms satisfies 

R<ma) e r 2 /p. (31) 

With use of (26), it follows that 

B<« f /^(«1/15). (32) 

If we neglect (his term we need consider only the simplified system 

p=-f m -ip (23.1a) 

pd= -ju e r 2 cos 2 $ sin 0. (23.2) 

9. Stage (2): Quasi-Steady Spiral Orbit 

It is convenient to treat stage (2) first. In this stage we shall assume that the change 4 in 
r in one revolution is so small that we can neglect the corresponding change in / compared to 
that in \f/. Presumably/ is proportional to the density, which in turn presumably varies like 
jR"exp[— (r— r e )/s] y where r e =mean earth radius=radius of sphere of same volume as the earth 
= 3,959 miles and where s is of the order 10 miles. Our neglect of the variation in /' in one 
revolution is thus equivalent to assuming that the corresponding decrease in r is somewhat less 
than one mile. (Clearly it would not do to make such an assumption for stage (3), when r 
may be diminishing much more rapidly.) 

To obtain the secular change, we can then "secularize" the right side of (23.2) by replacing 
cos 2 \p by Y 2 , thereby obtaining, 

pd=-ifco e r 2 sm 6. (33) 

Division of (33) by (23.1a) then gives 

dd/dp=%ma> e r 2 p- 2 sin 0. (34) 
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Note that the unknown drag function lias now dropped out from the theory — in any explicit 
form. In stage (2) 

v=7'\p, (35) 

closely enough, so that 

p=zmrv, (36) 

by (26) and (35). Equations (30) and (36) then yield 

p=7n(GM) 1/2 r 1/2 . (37) 

For convenience, we now define g e by 

^GM=g e r e 2 , (38) 

where r e =3959 miles, the above mean radius of the earth, and where g e is a corresponding 
effective value of the acceleration of gravity at sea level, free from centrifugal effects. Here 

g e = 32. 224 ft/sec 2 (39.1) 

(r e /&) ^=805.4 sec. (39-2) 

Then 

p=mr e g e ^r'K (40) 

From (34) and (40) it follows for stage (2) that 

esc 6 d6=0.25u e r e -hj e - lA r V2 dr (41) 



and thus that 



X! 2 csc e ^= ln Si=f w [<*w»- <^*i < 42 > 



where r x and d x are the radius and the inclination of the orbit at the beginning of stage (2) and 
r 2 and 6 2 are the corresponding quantities at the end of stage (2). 

An Example : Suppose that the initial altitude in stage (2) is 200 miles and the final altitude 
100 miles. Then 

r 1 = 3959 + 200 = 4159 miles (43.1) 

r 2 =3959 + 100=4059 miles. (43.2) 

We find 

(ra/r.)»-(r 1 /r.)»= -0.03860. (43.3) 

Then since co e =27r/86,164 sec -1 , it follows from (42), (43.3), and (39.2) that 

Thus 2 does not differ much from 6 X . If we now put 

<2^f (r e /fif e )^[(r 1 /r e )«-(r 2 /r e )^]>0 (45.1) 

d 2 /2 = v (45.2) 

ft/2— no, (45.3) 
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then, by (42) 

tan 17= exp (— Q) tan 170. (46) 

To find from (46) a power series in Q for 

A0 = 2 -ft (47) 

that will be valid for all ft, 6 2 in the interval (0, 7r), we use the fad that rj and rj are almost equal 
to each other, with each King in the first quadrant, and proceed as follows. Multiply (46) 

by i cos 7] and add cos 7j to each side. 
The result is 

cos rj sec 7] exp ir}=cos rj +i (sin r; ) exp (-Q) (48.1) 

= (exp ivo)[l —iy sin ^ exp (— i Vo )] } (48.2) 

where 

7^1-exp (-<?), (48.;]) 

a small positive quantity. Then 

17— i7 =Im In [1— iy sin %exp ( — ^ )] (49) 

—7 sin 77 cos 77 — 0.57- sin 2rj sin 2 ?7 +. . . , (50) 

on expansion of the logarithm and selection of the imaginary part. Thus 

A^2(r 7 -r 7o )=-7sin^,-7 2 sin ft sin 2 (ft/2) + . . . . (51) 

With use of (48.3), this becomes 

A6=-Q sin ft+Q 2 [0.5-sin 2 (ft/2)] sin ft + . . . , (52) 

which may be shortened, accurately enough, to 

Ad= -Q sin *i=-~ (re/ge) H [(ri/r e )K- (rjr § )*] sin ft. (53) 

Then, since 

^ (r 6 /^) H =0.009789 radians =0.5609°, (54) 

we have, as our final result for stage (2) : 

A0= -0.5609[(ri/r fl )*- (r,/r,)*] sin ft degrees, (55) 

where r € = 3959 miles. 

For the above example, where the altitude is assumed to diminish from 200 miles to 100 
miles in stage (2), insertion into (55) of the numerical value 0.03860 from (43.3 gives the result 

Ad=— 0.0217 sin ft degrees. (56) 

This is a very small effect. We next consider stage (1), where we can often expect a larger 
effect, but where the calculations are more uncertain. 

10. Stage (1): Shrinking Elliptic Orbit 

In stage (1) the density and thus the function/ in (23.1a) and (23.2) are appreciable only 
when the satellite is very close to perigee. In effect dd = and dp = unless r~r p . If we 
divide (23.2) by (23.1a), the result is 

esc 6 dd=mw e r 2 p~ 2 cos 2 \pdp. (57) 
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In stage (2) we denote initial quantites by the subscript and final quantities by the subscript 
1. The integral of (57) over the whole of stage (1) is then 

1 esc 6 dO=\n - ; 2, ; = mco e * r2 P~ 2 cos2 ^dv. (5$) 

>o tan (0 o /2) Jpo 

In the light of the above remarks we may put 

r 2 cos 2 \p=r 2 p cos 2 \[/ p (59) 

for any appreciable increment dp, with the result 

, tan (0,/2) pi 2 _ 2 

111 taiTW2) = mW Jp ^ C0S ^ (60) 

Now the perigee radius changes little even in stage (1), so that we shall not be far wrong if 
we replace r p in (60) by either r p0 or r\. Let us therefore replace r p by the arithmetic mean 

r=Ur PQ +n). (61) 

In (60) \f/ p is the angle between the line of nodes and the radius to perigee (the line of 
apsides) . When oblateness is taken into account (but not drag) , the rate of change of \p p is given 
by Garfinkel (see footnote 2, eq. 53), whose result in the present notation is 

^ p =3e(r e /l) 2 n(l-l.25 sin 2 9), (62) 

where 

e=(I p -I t )/Mr 2 e =0.00l09, (63) 

the quantities I p and I t being respectively the polar and transverse moments of inertia of 
the earth. We have here replaced Garfinkel 's equatorial radius by our mean radius r e , an 
unimportant change. The quantity / is the semi-la tus rectum, 6 is the inclination of the 
orbit, and n is the mean angular velocity of the satellite in its orbit; i. e., 

n=(vL/a 3 )X, (64.1) 

where the semi-major axis a is given by 

a=l{l-e 2 )- l =h(r a +r v ). (64.2) 

With the use of (64), (62), (38), and (39), we find 

^=36(^/r,)H(r»^(l-6 2 )- 2 (l-1.25sin 2 ^) (65.1) 

= 20.1°(r e /a) 7 ^(l-6 2 )- 2 (l-1.25 sin 2 6) per day. (65.2) 

As an example, suppose that r p0 = (3959 + 400) miles and r a0 = (3959 + 2000) miles and 
let us guess that r P i = r a i = ri= (3959 + 200) miles. (This latter figure corresponds to the 
beginning of stage (2), as assumed in our example above.) Then, by (64.2) and (29.2), a =5159 
miles and <? =0.1552, so that 6o=0.0241. Also ai=4159 miles and ei=0. The factor 
(r e /a) 7A (l — e 2 )~ 2 then increases from 0.415 to 0.842, with a mean value of about 0.628. Thus 

(^)mean-12.6°(l-1.25 sin 2 6) per day. (65.3) 

This amounts to 12.6° per day for an equatorial orbit or —3.15° per day for a polar orbit. 
Thus for a polar orbit the period of \p p is about 3.7 months. If the lifetime of a polar 
satellite is less than this, it remains doubtful what to do about cos 2 \l/ p in (60), although even, 
then we can establish an upper limit for |0i — 0o|. 
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If stage (1) lasts for a year or more, then \l/ p goes through a number of cycles in this state, 
unless sin 2 8 is too close to the value 0.8, for which the inclination 0=63.4°. If we consider 
various possible orbits and let approach 63.4°, the period of >fr v becomes greater and greater, 
until at this critical inclination the angle yf/ p would remain constant at its initial value during 
the whole of stage (1). 

For durations of stage (1) longer than a year and unless — 63.4°, however, we may find 
the secular change in d by replacing cos 2 \f/ p by ^ in (60), with the result 

ln !n^ U ^ .. . «*> 

= -Q.5o> e (rJgy>{rlr e )\trJn)y>- {rjk)*l, (67) 

with the use of (25) and (38). Here 

fr^O.SfoV+r-o 1 ) (67.1) 

from (29.1). 

We may evaluate d Y — O from the logarithmic tangent expression, in the same manner 
that we found 6 2 — 6\ in stage (2). From (67) we find 

0! — O =— 0.5Z, (0^63.4°) (68) 

and from (60) we find 

e l —6 = -L cosV„ (0-63.4°, with \p p effectively constant). (69) 

In (68) and (69) 

L - o> e (r e /g e ) »(f/r e ) 2 [{rjn) »- (rjl ) *| . (70) 

In any case 

|0 1 -0 O |<Z. (71) 

An Example: We consider the case where the initial perigee and apogee altitudes are respec- 
tively 400 miles and 2000 miles and assume that the final perigee and apogee altitudes of stage 
(1) are both 200 miles. The end of stage (1) in this example is then the same as the beginning 
of stage (2) in the previous example, so that the two examples together constitute one com- 
pound example. The r p0 =3959 + 400=4359 miles, ^=3959 + 2000 = 5959 miles, r pl = r al = 
ri = 3959 + 200=4159 miles, so that L=0.346°. Then, by (68) through (71), we have for 
Stage (1) (In miles, perigee 400-^200; apogee 2000 ->200) 

0i-0 o =— 0.173° sin O (6>o^63.4°) (72) 

0i — 0o=— 0.346 cos 2 ^psin O (0 O so close to 63.4° that \p p is effectively constant.) (73) 

In any case 

!0i-0o|<O.346° sin O . (74) 

These results may be compared with those for 
Stage (2) (In miles, altitude 200-^100) 

e 2 — l= — 0.022° sin ^ « — 0.022° sin O 
The over-all results are that 

02-0o= -0.195° sin O (0 O ^63.4°) 

|02— 0o|<O.368° sin O , in any case. 
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11. Summary 

The motion of an earth satellite relative to the earth's atmosphere produces a secular 
change in the inclination of the orbit. In predicting this change it is assumed that atmospheric 
drag acts in the direction of the air velocity relative to the satellite and that its magnitude 
depends on the relative air speed and in such a way on the altitude that it is appreciable only 
near the perigee, in the elliptic stage of motion. No other assumption is made about the drag. 

If d is the inclination of the orbit, p the orbital angular momentum, m the mass of the 
satellite, o) e the angular speed of rotation of the earth, r the distance of the satellite from the 
earth's center, and \p the angle between the radius vector and the line of nodes, it is then 
shown that esc d dd=mco e r 2 cos 2 \f/dp y to a good approximation, independently of the magnitude 
of the drag. 

The motion is separated into an initial elliptic stage, a quasi-steady spiral stage, and a 
final ballistic stage. The secular change in 6 is deduced for the spiral stage by integrating the 
above differential equation with use of an average value for r and the average value }{ for 
cos 2 \f/. For the elliptic stage the secular change is found by putting r 2 cos V equal to the perigee 
value r%cos 2 \p p in the above equation and then using some results of Garfinkel (see footnote 2) 
concerning the motion of the perigee. Nothing can be said about the final ballistic stage, 
however, without accurate knowledge of the drag. 

The secular change Ad in the inclination, a quantity difficult to measure, is thus expressed, 
for the first two stages of motion, in terms of easily measurable satellite distances. As an 
example, it is assumed that a satellite, with initial perigee and apogee altitudes of 400 miles 
and 2,000 miles, falls to an altitude of 200 miles at the end of the elliptic stage and an altitude 
of 100 miles at the end of the quasi-steady spiral stage. If these data should happen to be 
correct, then A0 would be —0.022° sin O for the spiral stage and —0.173° sin O for the elliptic 
stage. The latter figure depends on the assumption that 6 is not close to the critical value 
63.4°, for which there is no motion of the perigee relative to the line of nodes. In any case 
the total |A0| would always be less than 0.368°, for the assumed satellite altitudes, before the 
ballistic stage sets in. Once the ballistic stage begins, no predictions can be made by the 
methods of this paper. 

Washington, September 15, 1958. 
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